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^ . 1. Introduction 
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Thermal expansion and Hooke's law are among the fundamental properties of macroscopic 

matter. In many cases it has satisfactory microscopic explanation in the framework of equi- 
£f~) ■ librium statistical mechanics. This is a standard fact for the thermal expansion of gases and 

liquids. Derivation of linear theory of elasticity modules in general translation-invariant case 
see in (2 [3] . 

However, for crystals, where there is no translation invariance, the corresponding rigorous 
theory does not still exist. The nonstandard character of this problem is related, as we ex- 
plain below, to the coordination between local and global coordinate systems. The goal of 
this paper is to consider the simplest natural one-dimensional model in the global coordinate 
system. This model can be considered as one-dimensional crystal or one-dimensional network 
(macromolecules, microtubules etc.) in the biological cell, see 

Standard courses of physics, for example [BJ El E] ; or even more specialized books, for example 
|10l [TT] try to expalin thermal and elastic expansions of srystals via the oscillation of atoms 
about their equilibrium positions in the crystal lattice. In particular, the so called harmonic 
approximation and its nonlinear analogs are normally used. Any atom is considered to be in the 
vicinity of the equilibrium point. This is quite natural in the (microscopically) local coordinate 
system. However if we assume this for ANY atom in global coordinate system, then we cannot 
get any macroscopic expansion. 

The main idea is that the system is swelling uniformly at each point, and in any global coor- 
^ dinate system almost atoms become macroscopically far from their equilibrium positions. Note 

that similar situation was described in the so called random grammars, see |12] , Technically, 
this is achieved with another idea: to give atoms a possibility to be far away from their equi- 
librium position one should impose some restrictions on their motion. The simplest restriction 
would be to disallow particles jump through each other, that is the order of particles should be 
fixed. This idea of restricting the configuration space is close to one discussed by O. Penrose 
[1] in more complicated two and three dimensional cases. He argued that restrictions on the 
configuration space are necessary to get dependence of the free energy not only on the volume 
but also on its form. In our one-dimensional model this idea finds its rigorous justification 
together with very precise calculations. 

It is important to note that due to our restrictions one can give sense to the partition function 
of finite number of particles in the infinite volume. 

We consider the system of N + 1 particles (molecules, atoms) on the real line 

= x Q (t) < ... < x k {t) < ... < x N (t) 

where t > 0. The dynamics of this system is defined via smooth symmetric two-particle potential 
V(xi — Xj)' We make the following assumptions concerning this potential: 

(1) For u > the potential V(u) has a unique minimum at some point a > 0. Moreover, 
this minimum is quadratic; 

(2) We assume that 

V{0) = oo 

and make some comments on this assumption. In more realistic situation, when the 
real line is, for example, the x-axis in R d , the particles can deviate from this axis in 
perpendicular directions, and can in principle pass through each other in the x-direction. 
However time scale for such transitions is slower than their movement without changing 
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the order. Thus our assumption is a stable approximation to the realistic metastable 
situation. 

(3) The strongest technical assumption is that is we consider only nearest neighbor inter- 
action, that is the total energy is 

N 

v(x k - x k -i) 

fe=i 

This is quite natural for the one-dimensional networks but less natural for crystal mod- 
els. From this assumption it follows, in particular, that there exists a unique state with 
minimal energy (ground state) 

Xk = ka 

Without this assumption the situation with periodic ground states is more complicated, 
see 01(5], but we will not need this; 

(4) Our last assumption is that 

V(u) -> oo 

sufficiently fast (it is sufficient like > 1) as u —> oo' This assumption is necessary 

to make the partition function, for finite number of particles on the real line, finite. 
This is also a stable approximation to the metastable situation - to make impossible 
breaking the chain of particles into separate parts. 

2. Thermal expansion 
Assume that the left coordinate is fixed xq = and put 



Ui = Xi - Xi-i > 0, i = 1, N 
One can write the Gibbs density of the vector (tti, un) for temperature T = [i^ 1 as the 
following expression, which allows factorization, 

N N 

Z^ 1 exp(- fi^V (v,j))dx\...dxN = Z^ 1 exp(-/3V(ui))dui 

i=l i=l 

It follows that the random variables itj are independent, identically distributed, and do not 
depend on N. This makes the thermodynamic limit as N —¥ oo trivial. In particular the means 
< Ui > do not depend on TV, i and can be written as 

(2 1) <u >= m(T) - ^ uex P(~P V ( u ) du 

J exp(—/3V(u)du 

Then the total length of the chain is 

< x N >= Nm(T) 

We will consider here two cases: low temperatures, that is small deviations from the ground 
state and small perturbations around the state with fixed temperature. We will find conditions 
on the potential when macroscopic thermal expansion holds. 

2.0.1. Low temperatures. Assume that the potential has the following Taylor expansion at the 
point x = a 

V(a + y) = V(a) + c 2 y 2 + c 3 y 3 + c 4 y 4 + o(y 4 ),c 2 > 
Theorem 1. For low temperatures T = fi^ 1 

m(T) — a + m±T + o(T) 
where we will call mi the coefficient of thermal expansion, If C3 7^ 0, then 

3c 3 

This coefficient is positive iff C3 < 0. If on the contrary C3 = 0,C4 7^ 0, then 



" ll = 4c! 
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mi = 

Remark: If in the vicinity of a the potential is purely quadratic, that is it equals [x — a) 2 , 
then also m\ = 0. 

The physical example of the case C3 < can be the Lennard- Jones potential (except of the 
sufficiently large values of the argument) 

v LJ = C-^-C-f 

r r 

where (for a = 1) at the minimum a — \[2 we have 

63 

C3 = -V=2 

Proof. Consider the asymptotic expansion (for large /3) of the numerator and the denominator 
in (|2.ip . The standard Laplace method can be used. Then 

ex P (-(3V(u)du = J -exp(-l3V(a))(b + b 2 (2(3)- 1 + 0(/T 2 )) 

1 1 

uexp(-pV(u)du = J^exp(-pV(a))(d + d 2 (2(3)- 1 + 0(/T 2 )) 

The coefficients can be found from the equations 

00 

(2.2) x '( z) = ^6 fe z fc 

k=0 

or correspondingly 

00 

(2.3) (a + X {z)) X '(z) = Y. d * zk 

k=0 

where the function in the vicinity of the point z — 0, is defined by the conditions 

(2.4) V(a + x(z)) = V(a) + z 2 , x (0) = 0, x'(0) > 
see theorem 1.2.4 in [6j. Put 

X (z) = f lZ + f 2 z 2 + ... 

Then by flH|) 



2 

, c2.fl + 3c 3 f 2 f 2 + gff 

13 2/i 

Then from ([HI) , we get 

bo = fx,b 2 = 3/3 



That is why 



and 



The theorem follows. 



d a = af 1 ,d 2 = 3(/i/ 2 + af 3 ) 
dobo 1 = a 



_ d 2 d b 2 _ d 2 ab 2 _ 3 
TOl ~ Zb~Q ~ ~ZbJ ~ Zb~ ~ 2bo~ ~ 2 
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2.0.2. Small perturbations. One can write 

m(- ) = m(T = P^ 1 ) + m\e + o(e) 

P ^ 

and the differentiation gives 

mi =< uV >p — < u >/3< V >p 

It follows that the criterion of the positive expansion is the positivity of the correlation between 
random variables u and V. Again, if C3 7^ 0, then for large (3 the latter correltion is positive iff 
c 3 < 0. 

3. Hooke's law 

Hooke's law for the change Al of the length of cylindrical bar under the action of the force 
F is the following 

F Al 

~s ~ K ~T 

where S - the area of the cylinder crossection, I - is the length of the bar. 

Under the previous conditions assume also that there is an external potential V ex t (x) = —Fx, 
acting on the rightmost particle xn- Remind that xo = is assumed to be fixed. Then 

N N 

exp{-(3j2 V ( u r) - PV ext {x N )) = exp(-/3^(F( Ut ) - F^i)) 
i=i i=i 

as xn = ui + ••• + Define 

= / °° uexp(-(3(V(u) - Fu))du 
m[ ' ' f™exp{-/3(V{u)-Fu))du 
Then the elastic expansion is 

m(T,F) - m(T,0) 

Theorem 2. For fixed temperature T, sufficiently small F and any potential V, satisfying the 
conditions in the introduction, the elastic expansion is positive for positive F and negative for 
negative F. Moreover, 

m(T, F) - m(T, 0) = RF + o(F) 

where the elastic modulus 

(3.1) R - dm f F F) \ F = = f3(<ui>-<u > 2 ) 

Proof. The formula (|3.1[) is proved by simple differentiation. It is sufficient to observe now 
that the right hand side of (|3.1[) is the variance of nonzero random variable. 



4. Remark about the absence of expansion for oscillations around ground 

state equilibrium 

In the standard courses of physics normally the so called harmonic approximation is used 
with the hamiltonian (see for example [7], section 2.14) 

N N 1 

k,m=l k=l 

with symmetric matrix Cfc m , where = Xi—ia are the deviations from ground state coordinates. 
There is however an evident statement concerning joint distribution of the random variables . 
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Proposition 3. Let the distributions o/£i(t) are such that for any t, i and any sufficiently large 



P(Ut)\ >x) = oix" 1 ) 
Then for N — > oo and any t we have 

X(t) = maXiXi(t) — miriiXi{t) ~ Na 

In fact, for any e > the union of the events 

A i = {\t i (t)\>eN} 

has the probability which tends to zero as N — > oo. 

Thus, one cannot get any macroscopic expansion in this way. 
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